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ABSTRACT 


The main motivation for this thesis study is that 





Significant workload for aging transport aircrafts is 
related to dent removal from fuselages. This thesis is a 
preliminary investigation of aircraft fuselage dents using 


the Finite Element Method (FEM) via FEA ABAQUS software. We 








investigated single impact dent on fuselage panel at various 





locations and impact speeds. The material used for our 





finite element models is Aluminum Alloy 2024-T3, a typical 





material used for fuselages of older transport aircrafts. 


Our finite element model consisted of impact analysis, 





buckling prediction analysis, and postbuckling analysis 


successively. These analyses wer performed on both 








stiffened and unstiffened aluminum panels. We found that, 


depending on dent status in aluminum panel, the dent may 





increase or decrease buckling load of the panel compared to 


that of the virgin panel (undented). The buckling load of 





panel with low velocity impact is generally lower than that 





of the virgin plate. As the impact velocity is increased, 





buckling load of dented panel increases exceeding buckling 








load of virgin plate. In addition, we also noticed an 


existence of critical impact velocity at which the buckling 





load of the dented panel reached maximum and after which 





will start to decrease. 
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Ds INTRODUCTION 


A. BACKGROUND 


With today’s rapid increase in number of aging 








aircraft, the aerospace industry is fighting to reduce 
costs on maintenance and repair. Life extension of 
aircraft structures has become a major focus of the 


industry. The two main costly problems which have caught 





the attention of the scientific community are the fight 
against fatigue and corrosion [Ref. 3]. Much studies have 


been and are actively being conducted on crack fatigue, 





multiple site damages (MSD), aircraft skin inspection 








technology and amongst others [Ref. 4]. All these studies 





are, in one way or another, related to fatigue and 


corrosion crack damages of aircraft fuselage. Although a 





tremendous workload is also being spent on dent repairs and 


removal in large transport aircrafts [Ref. 1], not much 





study have been done on investigating the static stability 


of dent or impact damages on fuselage panels. 


After an extensive literature search and review, only 








a few studies have been found to involve impact damages on 
airplane structure. One in particular, by Cornelis Guijt 


from US Air Force Academy, performed experimental and some 





computer modeling study on fuselage dents. The study used 
a special impact swing hammer to develop various dent sizes 


on Al-2024-T3 fuselage panel. Although this paper focused 





mainly on fatigue investigation, it did present some static 





stability data for comparison. 


B. SCOPE OF RESEARCH 


This paper is a preliminary study on fuselage dents. 











The ultimate goal is to be able to provide guidelines on 
fuselage dent repairs. This paper studies the effect of 
fuselage dents on compressive failure loadings using 
computer simulation. The fuselage dents analyses were 


modeled using a computer simulation program called 
ABAQUS/CAE for both unstiffened and stiffened panels. After 


the deformed fuselages were created, they were subjected to 





compressive loading and eigenvalue calculation analyses. 
The buckling load results between the compressive failure 
analyses and eigenvalue analyses were compared. The 
buckling behaviors of unstiffened and stiffened panels were 


also presented. 


II. MODELING 


A. UNSTIFFENED FUSELAGE PANEL 


The unstiffened fuselage panel was modeled using 


ABAQUS/CAE as a 0.508 m by 0.508 m (20 in. by 20 in.) 3-D 





deformable shell panel. The shell panel was created via 
extrusion of one of the x-axis edge. The material property 


was defined as elastic-plastic with shell thickness of 








0.0015875 m (0.625 in.). Density of the shell plate was 
2780 kg/m3, elastic modulus was 73.1 GPa with Poisson ratio 


of 0.33. An 80x80 S4R (quadrilateral shell) finite element 





mesh was used for the shell plate. A summary of the 


material property and geometric dimension are tabulated in 








Table 1&2 and a finit lement mesh of the fuselage panel is 


given in Figure l. 




































































Table 1. Material Property of Shell Plate A1-2024T3 
(Ref. 5) 
Young's Modulus of Elasticity T34L GPA 
Yield Strength 3.45E+08 Pa 
Ultimate Strength 4.83E+08 Pa 
Mass Density 2780 Kg/m?> 
Elongation at Failure 15% 
Shell Plate Dimension (LxWxt) 0.508x0.508x0.003175 mxmxm 
Table 2. Geometric Dimensions of Fuselage Panel (after 
Ref. 1) 
Model Geometry 
Unstiffened Plate Dimension 0.508x0.508 m 
Shell Thickness 0.15875 cm 
Stiffened Plate 
Dimension 0.508x0.508 m 
Shell Thickness 0.0015875 cm 
Stiffener Dimension 1.905x1.905x1.905 cm 
Shell Thickness 0.15875 cm 
































Figure l. Finite Element Mesh of Shell Panel-(element 
size=0.635 cm). 


B. STIFFENED FUSELAGE PANEL 


The Stiffened fuselage panel was modeled using 
ABAQUS/CAE similar to the unstiffened panel. The main 
difference was the addition of hat shaped stiffener 
supports. The 0.508 m by 0.508 m (20 in. by 20 in.) 
stiffened panel consisted of a set of hat stiffeners at 
20.32 cm (8 in.) apart. Fach of the joint between 
stiffener and main panel was modeled to be twice as thick 
as (0.3175 cm or 0.1250 in.) the main panel itself. The 
mesh size was the same as the unstiffened panel but 
included the mesh for the two stiffeners. The height and 
width of winged portion of each stiffener were both 1.905 
em “(Og 75: ins) < The material property used was same as 
listed in Table 1. A detail dimension of stiffened panel 


is shown in Table 2. 





Figure 2. Finite Element Mesh for Stiffened Panel 
(element size=0.635 cm). 


Cc. UNSTIFFENED FUSELAGE PANEL WITH PROTRUDING DEFORMATION 


An additional study was conducted to investigate the 
stiffening effect caused by the present of fuselage dents on 
compressive failure loads. The model was created using 
ABAQUS/CAE similar to the unstiffened fuselage panel. The 
model consisted of a 0.508 m x 0.508 m deformable shell 





panel with protruding hollow cone-like structure as shown in 





Figure 3. The deformable shell panel with 20.32 cm (8 in.) 





diameter center cutout was created separately from the 





protruding structure. ABAQUS has a feature in which parts 
with same uniform material property can be merged or 
combined into one piece. In this case, the cutout shell 


panel was merged with the protruding structure to form our 





final modeling part. The material property used for the 


protruding structure was same as the rest of the shell panel 
(Table 1). Figure 4 showed a side view of a protruding 


structure with the dimensional property. 





Figure 3. Finite Element Mesh of Shell Panel with 
protruding structure with 20.32 cm (8 in.) 
diameter cutout. 





Side view of protruding structure. 
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III. SIMULATION 


A. INTRODUCTION TO ABAQUS/CAE 


All computer simulations in this study were conducted 








using a finite element simulation program called ABAQUS. 
ABAQUS is an advanced finite element analysis that provides 


complete and powerful solutions for linear and nonlinear 








engineering problems. It is a suite of finite element 
analysis modules consisting of preprocessor, solver, and 
postprocessors. Preprocessor allows the creation and 
assembly of modeling parts. Their solver package includes 
ABAQUS/Explicit and ABAQUS/Standard. ABAQUS Viewer is the 
postprocessor which allows visualization and post processing 


of analysis results. Their most complete package, 





ABAQUS/CAE, is a fully interactive environment equipped with 














graphical user interface (GUI) of menus, icons, and dialog 
boxes. It provides the complete working environment for 
users to create models, submit jobs for analysis via the 





solver package, and view and process results [Ref. 2]. 


B. MODELING AND SIMULATION PROCEDURE 


Both the stiffened and unstiffened plate models were 
subjected to three different analyses. The first analysis 


was the dynamic impact simulation in which an artificial 





dent is created via a 5.08 cm (2 in.) diameter rigid ball at 











various locations and impact velocity. The deformed plates 





were than imported along with the material states 


(stress/strains) using ABAQUS/CAE in order to undergo 





compressive failure analysis, which was conducted in two 


different ways. The first study was the incremental 





compressive loading to the impacted panels until the 








buckling of the panels. In this study, the residual stress 





and strain caused by the impact was included in the model as 


well as the dent shape. The other study was undertaken for 





eigenvalue analysis of the dented panels. In the latter 








case, the deformed shapes of the panels were considered in 


the analysis, but not the residual stress and strain. Both 





results were compared to the buckling loads of the virgin 





panels before impact. 














A total of 21 impact simulation runs were conducted and 
each took average of 15 hours. Fach of the 28 compression 
test simulations took from 1-3 hours to complete. Finally, 











the time required for each eigenvalue calculation analysis 





was approximately 30 minutes. All analysis were conducted 











using a Dell Dimension Desktop equipped with single Pentium 
4 processor at 2.4 GHz and 1 GB RAM. A summary of all 


simulations is given in Table 3. 





Cc. DYNAMIC CONTACT IMPACT SIMULATION 





The ABAQUES/Explicit analysis engine was used to 
perform the dynamic impact simulation for each fuselage 
model. ABAQUS/Explicit possessed a powerful nonlinear 


dynamic analysis capability which allowed us to define 




















contact interaction. In this analysis, surface-surface 
contact interaction was defined for the rigid ball impactor 
and the deformable fuselage panel. The normal behavior of 





interaction property was defined as hard contact. Contact 
pair separation after impact was allowed in order to 
eliminate contact load after impact on the deformed fuselage 


panel. The frictionless formulation was used for the 
10 





tangential contact interaction behavior. The Kinematics 





Contact Method was used as the mechanical constraint 


formulation. The boundary conditions of the fuselage panel 











were rigidly fixed on all four edges with an initial 











condition of the model defined for the rigid ball impactor 
velocity. An enhanced hourglass mesh control was chosen in 
order to reduce th ffect of hourglass stiffness. This 





mesh control property propagated into subsequent analyses 
and did not require further adjustment. For each fuselage 


panel model, a series of impactor velocity was used to 





investigate the various dent sizes. Total simulation time 





step was set at 1 second. A complete detail of the creation 





of this analysis model is provided in Appendix A. 
D. BUCKLING PREDICTION/EIGENVALUE CALCULATIONS 


ABAQUS/Standard analysis engine was used for buckling 








prediction analysis. Individual deformed geometry created 





from the impact analyses was imported into a new model to 





obtain eigenvalues for the first 5 modes. A simple support 





boundary condition was used on the two opposite ends while 











the other two opposite ends remained free. The in-plan 


compressive loading was applied to one of the simply 


supported edges. A uniform shell edge load of 1 Newton was 
used. The total predicted buckling load is found by 
multiplying th respectiv igenvalue value by 1 Newton. 





The eigensolver selected for this analysis was Subspace 








iteration method. A detail description of the creation this 


model is provided in Appendix A. 
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E. COMPRESSION FAILURE ANALYSIS/POSTBUCKLING SIMULATION 


ABAQUS/Standard analysis was also used for compression 


failure analysis. Unlik th igenvalue calculation 














analysis, th deformed fuselag geometry along with 














respective residual stresses and strains were imported into 
a new model in order to perform the compression test 
simulation. Boundary conditions used were simple supports 
as described in the previous’ section. An incremental 


uniform shell edge load was used for all compression test 





analyses. Static Riks Iteration Method was chosen for 





buckling calculation. 


Table 3. Fuselage Panel Simulation Summary. 





Al-2024-T3 Fuselage Panel Simulations 


Impact Type 
Panel Type 


ndeformed Panel 


ndeformed Stiffened Panel 








Deformed Unstiffened Paneljpal 


Center 





Center 
Center 





Center 





Center 





Center 
Center 





Deformed Stiffened Panel Center 





Center 





Center 





Center 





Center 





Center 





Center 





stif 





stif 





stif 





stif 





stif 





stif 








L L L L L L L L L L L L L L i L L L L L 
t t t t t t t t t t t t t t t t t t t t 
Pree RL eee elelelelelelelelelelelebelre 











stif 
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IV. RESULTS/DISCUSSIONS 


A. DYNAMIC IMPACT SIMULATION 


The impactor chosen for dynamic impact simulation was a 
0.0254 cm (1 in.) radius rigid ball. The rigid ball 
impacted a shell panel at various velocities and locations 
(Table 3) creating different dent sizes and depths. As 


expected, the dent sizes and depths increased with the 





increasing impact velocity. Dent sizes varied from 6.35 cm 
(2.5 in.) to slightly less than 15.24 cm (6.0 in.) in 
diameter as shown in Table 4. The depth of dents varied 
from about a 0.3175 cm (1/8 in.) to a little over 6.35 cm 
(2.5 in.). Although dent sizes and depths were not the same 
for the deformed stiffened and unstiffened panels subjected 


to the same impact condition, both did show the same trend 








of an increase in dent sizes and depths with respect to an 


increase in impact velocity. At lower impact velocities, 








however, it was difficult to determin th dent radius 


because the dent sizes were too small to measure accurately 





with finit lement mesh size of 0.635 cm (0.25 in). For 








instance, with stiffened panel impacted at the velocity of 








10 m/s at the center, the dent size was probably less than 
0.635 cm and therefore too small to be measured. The 
existence of a dent was confirmed by the contour plot of 


equivalent plastic strain. 


ales: 








Table 4. Fuselage Dent Results 


Fuselage Simulation Dent Result Summary 


Impact Dent Dent 
Speed | Radius |Radius| p epth Depth 
(m/s) (em) (in) (in) 


Impact 
Location 
Panel Type 


ndeformed Panel 


ndeformed Stiffened Panel N/A 
Deformed Unstiffened Panel Center 








Center 





Center 





Center 





Center 





Center 





Center 
Deformed Stiffened Panel Center 





Center 





Center 


Center 





Center 





Center 





Center 


stiffener 





stiffener 





stiffener 





stiffener 





stiffener 





0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
O). 
0. 
O. 
0. 
O.. 
O. 
0. 
0. 
0. 
0. 
0. 


stiffener 





Figure 5 shows the plot of dent radius versus impact 
velocity. The plot showed that initially there was an 


increase in the dent size with an increasing impact velocity 





for all three cases. For the center impact of stiffened and 
unstiffened panels, the dent sizes showed a large increase 


from a low velocity (10 m/s) to a higher velocity (30 m/s) 





and slowed down or stabilized at a much higher velocity. The 





dent size of the stiffened panel with an impact at 10 m/s 





was too small to measur ven from the contour plot of 


equivalent plastic strain. 
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For the direct stiffener site impact, however, the 
figure showed the dent radius experienced a large jump 
initially and decreased at high velocities. There is no 
clear explanation for that at this time. This phenomenon 


needs to be further investigated. 





Dent Radius vs Impact Speed 








50 —®—Unstiffened Panel 





—@ Stiff Panel Center 
Impact 
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(ee) 
Oo 
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Figure 5. Dent size impact results. 
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Dent Depth vs Impact Speed 
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Figure 6. Dent depths impact results. 








The dent depths were determined by the maximum 
transverse displacement shown on the displacement contour 
plot of each deformed panel. Figure 6 showed that the dent 
depth increases monotonically as a function of impact 
velocity. Both the unstiffened panel and the direct 


stiffener impact panel showed a nearly linear relationship 





between the dent depth and the impact speed. However, the 
stiffened panel with the center impact indicated a less 


linear relationship. 


B. EIGENVALUE ANALYSES OF VIRGIN PANELS 


EFigenvalue analyses were performed for all virgin 








panels before impact in order to compute their initial 





buckling loads. The first five buckling modes and their 


associated buckling loads were calculated. However, the 





lowest buckling load was selected as the reference value to 
be compared with the compressive failure load after the 
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impact damage. The critical buckling loads for the 


unstiffened and stiffened virgin panels were 990 N and 





15,920 N, respectively. Some of the buckling mode shapes for 


unstiffened panels are shown in (Figure 7-10). 


1000e+00 





: Buckle _Plate.ocdb ABAQUS/STANDARD Version 6.6-1 Wed Mar 28 16:13:01 : 
3 Step: Step-1, Buckling prediction of undeformed plate. 
Mode 1: BigenValue = 989.65 
Primary Var: U, Magnitude 
Deformed Var: U Deformation Scale Factor: +5.080e-02 
Figure 7. Virgin unstiffened panel buckling Mode 1. 





















U, Magnitude 
+1.000e+00 
+9.167e-01 
+8. 333e-01 
+7.500e-01 
+6.667e-0O1 
+5.833e-01 
+5.000e-01 
+4.167e-01 
+3.333e-01 
+2.500e-01 
+1.667e-0O1 
+8.333e-02 

-000e+00 





3 ODB: Buckle_Plate.odb RBAQUS/STANDARD Version 6.6-1 Wed Mar 28 16:13:01 
Step: Step-1, Buckling prediction of undeformed plate. 
Mode 2: EigenValue = 2759.4 
Primary Var: U, Magnitude 
Deformed Var: U Deformation Scale Factor: +5.080e-02 
Figure 8. Virgin unstiffened panel buckling Mode 2. 


Lf 


U, Magnitude 
+1.000e+00 


*5.833e-01 
+5.000e-01 
+4.167e-01 





3 ODB: Buckle _Plate.odb ABAQUS/STANDARD Version 6.6-1 Wed Mar 28 16:13:01 
Step: Step-1, Buckling prediction of undeformed plate. 
Mode 3: EigenValue = 4069.5 
Primary Var: U, Magnitude 
Deformed Var: U Deformation Scale Factor: +5.080e-02 
Figure 9. Virgin unstiffened panel buckling Mode 3. 


U, Magnitude 
+1.000e+00 
+9.167e-01 
+8.333e-61 
+7.500e-01 
+6.667e-01 
+5.833e-01 
+5.000e-01 


+0.000e+00 





ODB: Buckle_Plate.odb ABAQUS/STANDARD Version 6.6-1 Wed Mar 28 16:13:01 


Step: Step-1, Buckling prediction of undeformed plate. 
Mode 4: EigenValue = 5908.5 

Primary Var: U, Magnitude 

Deformed Var: U Deformation Scale Factor: +5.080e-02 


Figure 10. Virgin unstiffened panel buckling Mode 4. 
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Cc. COMPRESSION TEST/POSTBUCKLING SIMULATION 


Following the dynamic impact analyses, the deformed 








shell panels including the residual stresses and strains 
were imported into a new model and were subjected to in- 


plane compressive loading. The compressive failure loads 





from the compression analyses were determined from the load 








and displacement data. For the unstiffened panels, critical 
failure loads can be read directly from their load versus 
displacement curves. For instance, Figure 11 showed load 


versus displacement curve for a deformed plate formed by a 


35° m/s Anpact. The curve clearly showed the critical load 
at 750 N. The failure modes for unstiffened panels were as 
expected (Figure 12). Their shapes were a half sine curve. 


The failure loads for stiffened panels, on the other 





hand, were not as obvious by looking at the load 











displacement curve (Figure 13). In this case, the failure 


loads were determined from the last data point in the graph 





where the finite element solution diverges. 
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Load vs Displacement (N 41) 
(Deformed Plate @ V=35m/s) 
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Figure 11. Load displacement for deformed plate at 
V=35m/s. 










8, Mises 
SNEG, (fraction = -1.0) 
(Avg: 75%) 


+9.578e+06 









2 Postbuckling analysis of defo: 


ODB: POSTBUCKLE_Plate_V35.odb 
3 Step: Step-1 


Increment 42: Arc Length = 533 
Primary Var: S, Mises 
Deformed Var: U Deformation Scale Fact 


=35m/s using ced h 
6.6-1 Mon MamO5 21:27:06 Pacific 





Figure 12. Shell panel after impact (35m/s) and 
compression tests. 


20 





Load vs Displacement (N 536) 
(Stiffened Plate Center Impact @ V=54m/s) 
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Figure 13. Load displacement for 
plate at V=55 m/s. 


deformed stiffened 










Ss, Mises 
SNEG, (fraction = -1.0) 
(Avg: 75%) 

+4.812e+08 





Postbuckling analysis of er) impact at v=35m/s: 

ee ODB: POSTBUCKLE_St iffenedPlate™ WARD Version 6.6-1 

3 1 Step: Step-1 
Increment 100: Are Length = 7.487 
Primary Var: 8, Mises 
Deformed Var: U Deformation Scale Factor: +1.000e+00 


Figure 14. Shell panel after 
compression tests. 


en 
Tue Mar 06 12:0% 


impact (55m/s) and 





The compressive test analyses for the damaged 
unstiffened panels showed that certain fuselage dents seemed 


to increase the failure load of the panel compared to the 
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buckling load of the virgin panel. An experimental study by 
Guijt also noted a strengthening effect caused by existence 


of dents (Ref. 1). From the summary listed in Table 4, it 








can be seen that at low impact velocities, the compressive 


failure loads were smaller than the virgin buckling load. 





In fact, at the impact speed of 10m/s, 30m/s, and 35m/s, the 





compressive failure loads were 606 N, 659 N, and 725 N, 





respectively, which were less than that of the virgin panel 


(990 N). As the impact velocity was increased, the failure 





load continued to increas xceeding the buckling load of 





the virgin panel. The failure load increased to a maximum 





at 1368 N and began to decrease. A plot of failure load 





versus impact speed is given in Figure 13. The failure 


loads have been normalized against the buckling load of the 





virgin panel. The plot showed a general stiffening effect 


caused by impacts at higher velocities. An explanation of 





this stiffening effect is given in the last section of this 





chapter. 


Table 5. Unstiffened Panel Buckling Load. 








Impact | Dent Dent | Dent | Dent 

Speed |Radius/Radius|Depth |Depth 
(m/s) (cm) (in) | (em) | (in) 
Center 10 3.18 1.25 |0.30]0.12 606 


Failure 
Load (N) 


Impact 
Location 

















Center 30 6.67 2.63 |1.52]0.60 659 
Center 35 7.58 2.98 |2.00]0.79 725 
Center 55 6.61 2.60 |3.14]1.24 1093 
Center 60 6.55 2.58 | 3.75 ]1.48 1206 
Center 65 7.00 2.76 |4.18]1.65 1368 
Center 70 6.75 2.66 |4.69)1.85 1351 
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Buckling Load at Various Impact Speed 
(Unstiffened Panel) 

2 1 
AS} 1.40 | 
2 1.20 | 
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iN 
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Impact Speed (m/s) 
Figure 15. Buckling load for deformed unstiffened panel. 


For the deformed stiffened shell panels, results from 





compression test analyses revealed similar behavior as 
earlier. The failure loads for damaged stiffened panels 
started out at small values and increased, surpassing the 
buckling load for the virgin panel (15,920 N), to a maximum 
and began to decreased back down. From Table 7, for the 


stiffened panel with center impact at 10 m/s, the failure 





load was 15,190 N. This was slightly lower than buckling 


load for the virgin panel. The failure load reached a 





maximum value of 42,260 N at 30 m/s and started to decrease 
rapidly to 15967 N at 70 m/s. These results, again, 
appeared to suggest a stiffening effect caused by high 


velocity impacts. 


For stiffened panels with direct damage to_ the 
stiffener site, it was found that the maximum buckling load 
was 39,230 N at the impact velocity of 55 m/s. At low 


impact velocities, the buckling loads were approximately 25 





%& lower than that of the virgin panel without dent. It 
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appeared that damage to the stiffener support at low 
velocity impact can be quite detrimental to the structure of 
a fuselage. The failure load versus impact speed plot in 
Figure 16 and 17 below showed a summary of buckling behavior 


of deformed stiffened panels. 














































































































Table 6. Stiffened Panel Buckling Load Results 
Summary. 
Impact Dent Dent Dent Dent ; 
Bae Speed Radius Radius Depth Depth earn 
(m/s) (cm) (in) (cm) (in) 
Center 10 0.27 0.11 15186 
Center 30 5.06 2.39 1,15 0.45 42261 
Center 35 6.07 2.70 1.41 0.56 39295 
Center 55 6.87 2.69 1.89 0.74 18044 
Center 60 6.83 2.75 2.23 0.88 16278 
Center 65 7.00 2.81 269 1.06 17690 
Center 70 cares) 0.00 iS eae ae 1.25 15967 
On stiffener 30 0.41 0.16 11459 
On stiffener 35 3.94 2.48 0.45 0.18 12084 
On stiffener 55 62.31 2.23 2.48 0.97 39232 
On stiffener 60 5.66 2225 2.86 ls LZ 32719 
On stiffener 65 5.72 1.66 3.16 1.24 28121 
On stiffener 70 4.21 0.00 3.47 1.36 28719 
Buckling Load at Various Impact Speeds 
(Stiffened Panel Center Impact) 
73 3.00 
@ 
° 
2.50 
D 
£ 
Ss 2.00 
(3) 
1.50 
3 
& 1.00 
£ 0.50 
-— 
° 
= 0.00 
0 10 20 30 40 50 60 70 80 
Impact Speed (m/s) 
Figure 16. Buckling loads for stiffened panel center 
impact. 
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Buckling Load at Various Impact Speeds 
(Stiffened Panel Stiffener Impact) 
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Figure 17. Buckling loads for stiffened panel stiffener 


impact. 


D. EIGENVALUE ANALYSES OF DAMAGED PANELS 





Linear eigenvalue analyses wer performed for each 
damaged fuselage panel after impact in order to compute 


their buckling loads. The geometric models of the damage 





panels were imported from the impact analysis models for 








eigenvalue analysis. Because the linear eigenvalue analysis 
does not require residual stress/strain, this information 
was not imported. The buckling loads of the damaged panels 
were compared with the compressive failure loads after the 
impact damage. Even though the boundary conditions and the 


loading direction were the same for the two analyses (i.e. 











nonlinear compressive and linear eigenvalue analyses), the 


comparisons showed a large difference in the failure loads 





and mode shapes between the two analyses as seen in Table 7. 





For instance, at 70 m/s impact on stiffener, the damaged 
stiffened panel showed a buckling load of 15600 N while the 


failure load from compressive test analysis gave slightly 
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over 28700 N. Moreover, when comparing buckling mode shapes 
for the same damaged panel, the compressive test analysis 
showed a larger deformation along the stiffener support 
(Figure 18). On the other hand, the linear eigenvalue 
analysis showed a larger deformation on one of the free 
edges (Figure 19). The next four eigenvalue mode shapes 


for the same panel were also different from the compressive 





test analysis. The results for the compressive failure 
analysis and for the first mode eigenvalue calculation are 


shown in Figure 18 and 19 respectively. 



















8, Mises 
SNEG, (fraction = -1.0) 
(Avg: 75%) 
+5.687e+08 
+5.222e+08 
+4.757et08 
+4.292e+08 
+3.826e+08 
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+2.431e+08 
+1.965e+08 
+1.500e+08 
+1.035e+08 
+5.695e+07 
-042e+07 










Postbuckling analysis of defo ffener) at impact speed 
ple ODB: POSTBUCKLE_StiffenedPlate_O RNDARD Version 6.6-3 Sun Mar 04 13 


Step: Step-1 
1 Increment 100: Arc Length = 13.5 

Primary Var: S, Mises 

Deformed Var: U Deformation Scale Factor: 


Figure 18. Buckled stiffened panel with stiffener impact 
at 70m/s (Compressive Test Analysis). 
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+8.337e-01 
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+3.335e-01 
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Buckling prediction of defo: fer) impact V=70m/s; e 


ODB: BUCKLE_StiffenedPlate_Os_V Version 6.6-3 Sun Mar 04 14:09:26 F 
Step: Step-1 
3 Mode 1: EigenValue = 15611. 


Primary Var: U, Magnitude 
Deformed Var: U Deformation Scale Factor: +5.113e-02 


Figure 19. Results from Eigenvalue Calculation Mode 1. 


Table 7. Compare Eigenvalues and Buckling Loads. 


Fuselage Simulation Eigenvalues Result Summary 


ndeformed Panel 

ndeformed Stiffened Panel N/A 15922 
10 2125 
65 TLS 
10 8978 
25 20094 
55 20687 
60 23357 
65 18858 
70 6245 
70 15611 














Deformed Stiffened Panel 
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E. EXPLANATION FOR STIFFENING EFFECT OF FUSELAGE DENTS 


From the compressive test analyses, it was noticed that 





the existence of a dent seemed to cause a strengthening of 
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the deformed fuselage panels as discussed in previous 





sections. It was noticed that at lower impact velocity, the 


failure load was smaller than that of the undented virgin 








panel. At a higher velocity, the failure load rises above 





that of the virgin panel and reached a maximum at a critical 





impact velocity. At an impact velocity greater than the 








critical velocity, the buckling load began to decrease. 
This important observation was made in both the unstiffened 


and stiffened panels. 


By examining the dent sizes and depths, the dent shape 





was not significant at a low impact velocity. The dent was 
a very small localized deformation at the impact site 
compared to the rest of the panel. However, the local 
residual stress at the dent caused by the impact reduced the 


overall strength of the panel. As a result, a very small 








local dent reduced the compressive failure load lower than 





the buckling load of virgin plate without a dent. As the 


impact velocity increased, the dent shape became more 





Significant. Here, there were two competing factors. First 
the dented site has a lower local strength with high 


residual stresses which would reduce the compressive failure 





load. On the other hand, the shape of the dent played a 





local bending stiffening effect which would increase the 
compressive failure strength. A reduction in compressive 
failure load can be viewed as a negative effect while an 
increase in compressive failure loads a positive effect. 


Thus, up to a certain dent size and depth, the net effect 





(i.e. negative effect subtracted from the positive effect) 
increased positively. As a result, the compressive failure 


load of a dented panel increased with the impact velocity. 








The compressive failure load eventually exceeded that of the 
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virgin panel. After reaching a peak compressive failure 


load, the net effect would start to decrease resulting in 











decrease of the compressive failure load. 





In an attempt to support the above explanation, a 








series of finite element models were created. The fuselage 
panel with protruding structure model was described in the 
Chapter 2, Sec. C. This model was an attempt to investigate 


the behavior of a dented fuselage panel via eigenvalue 





analysis. The height of protruding structure modeled th 
dent depth while the siz of the cutout modeled the dent 








size. A series of eigenvalue analysis were performed using 


various protruding heights. Table 8 summarized the results. 





In the model, a cutout hole represents the reduced 





strength of the dented section caused by residual stresses 


while the protrusion represents the increase of bending 








stiffness resulting from a dent. As seen in the table, when 
the effect of the hole is greater than the protrusion 


effect, the buckling load is lower than that of the panel 





without a hole. However, as the protrusion effect (i.e. 


increase of bending stiffness with dent) becomes more 








dominant, the buckling load becomes much greater than that 


of the plate without a hole. 
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Table 8. Buckling Loads for panel with protrusion. 





Protrusion 
Height h 
(cm) 


Protrusion| Buckling 


(N) 


Panel without cutout 

8" cutout w/o protrusion 
8" cutout w/ protrusion 
8" cutout w/ protrusion 














8" cutout w/ protrusion 
8" cutout w/ protrusion 
8" cutout w/ protrusion 
8" cutout w/ protrusion 
8" cutout w/ protrusion 
8" cutout w/ protrusion 
8" cutout w/ protrusion 
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Vv. 


CONCLUSION/RECOMMENDATION 














This thesis investigated th ffect of fuselage dents 
on compressive failure load using computer modeling and 
simulations. By modeling the different impact velocities 
and locations, various dent shapes (sizes and depths) were 
created for both unstiffened and stiffened panels. An 
additional study was also performed to attempt to further 
explain the stiffening effect of fuselage dents on the 
compressive failure load. From this study, a few 
conclusions and generalizations were made based on 


examination of the results. 


Ts 





Depending on the dent status in a panel caused by 
an impact, the dent may decrease or increase the 
failure load of the panel compared to the virgin 
panel without impact. 








It was observed that at a low impact velocity, 
failure loads of damaged panels were generally 
smaller than that of the virgin panels. However, 
as the impact velocity was increased, the failure 
load also increased and eventually surpassed the 
buckling load of the virgin panels. At a certain 
critical impact velocity, the failure load reached 
a maximum after which it began to decrease. In 
general, th xistence of a dent can strengthen a 
deformed panel in a certain dent size range or 
impact velocity. 




















A direct impact or damage 
a low impact velocity can 
failure load was shown to 
of the virgin panel withou 


Th linear 





to the stiffener site aft 


be detrimental as the 
be much less than that 
t impact. 








igenvalue calculations of the damaged 


panels did not compared well to the compressive 


failure load of the sam 





stresses wer not 
analysis. 
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accounted 





residual 
eigenvalue 


panel becaus 
in the 





A few recommendations can be made for future studies 


based on the results obtained in this investigation: 


Te 








In reality, dents can be of many geometric shapes 
and sizes. Recommend using different impact 
shapes (i.e. a cylinder impactor at an angle). 


Only two impact locations were investigated in 
this study. Recommend investigating other impact 
locations to fully understand the effect of dents 
on stiffened panel. 








A hat-like stiffener model was used in this study. 
Based on literature reviews, other stiffener 
designs are available for fuselage panel. 
Recommend using a different stiffener design. 
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APPENDIX A.. STIFFENED AND UNSTIFFENED PANEL IMPACT 
MODELING USING ABAQUS/CAE 


This appendix explains in detail how impact modeling 


was created using ABAQUS/CAE. The final result would be an 





input file which is used to submit into a job request for 
analysis using ABAQUS/Explicit. ABAQUS is an advanced 


finite element analysis that provides complete and powerful 











solutions for linear and nonlinear engineering problems. It 





is a suite of finite element analysis modules. ABAQUS/CAE, 





having a modern graphical user interface (GUI) of menus, 





icons, and dialog boxes, provides the most complete 


interface with the ABAQUS solver programs available [Ref 2]. 





In order to understand the process, the user should be 
familiar with the components of the ABAQUS/CAE and the 


appearance of the window. Figure 20 shows the components 





that appear in the main window. 
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Title bar Menu bar Toolbar Comtext bar 
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area drawing area area command line interlace 
Figure 20. Components of the ABAQUS/CAE main window. 


The part (stiffened plate) for the model can be 
created in the part module of the main window. 


From the main menu bar: 


Part > Create > Name the part “Stiffened Plate” 
> Choose 3-D modeling space, deformable type, 
shell shape feature, and extrusion type feature > 


Set approximate size to 1.0 > Click OK (Figure 
21) 
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@® Create Part x! 
Name: [stiffened_Plate 


Modeling Space 















@ 3D © 2DPlanar © Axisymmetric 


Type Options 


© Discrete rigid None available 


© Analytical rigid 


Base Feature 









Shape Type 
™ Solid 
Shell 

Revolution 
Wire 

Sweep 
™ Point 


Approximate size: fi 
Cancel _| 


Figure 21. Create Part Dialog Box 








Dick From the edit tools to the left of the viewport, 





select the “+” (create isolated point tool) to 
enter following coordinates: (-0.254,0) and 
(0.254,0). Select “create lines tool” (to the 
left of create isolated point tool) and connect 


the two points previously created. Follow the 





specification given in Table 2 to create two hat 
stiffeners at 0.2032 m (8 inches) apart using the 
“create isolated point” and “create lines” tools. 


Once complete (see Figure 22): 


Click Done at the bottom of viewport > in Edit 
Base Extrusion box, enter 0.508 for depth > Click 
OK (Figure 23) 
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Figure 22. Line part shell extrusion. 
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ie 
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Figure 23. Unmeshed Stiffened Panel 
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Figure 24. 


The rigid ball impact is created in a similar 


fashion as earlier. From the main menu bar: 


Part > Create > 3D Modeling Space > Analytic 
Rigid Type > Revolved shell > Approximate size 1 
> Click Continue. > Select Create Arc tool > 
Follow direction at bottom of viewport; select 
center point at origin; select perimeter point at 
(0,-0.0254) and second point at (0.0254,0). 73 
Create second arc with at origin and perimeter 
points at (0.0254,0) and (0,0.0254). > Click OK. 


Create reference point at top of ball > Tools > 
Reference Point > Select top point of ball 73> 
Click Done. See Figure 24. 








Rigid Impact Ball with reference point. 
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The 


material property is assigned to the stiffened 


plate in the Property Module. First, the material 


property must be defined. From the main menu bar: 


Figure 25. 


Material > Create > Name material “A1-T3” > In 
General pull down menu, select density, enter 2780 
kg/m? for mass density. > In Mechanical pull down 
menu, select Elasticity > Elastic > Enter 
Young’s Modulus and Poisson’s Ratio provided in 
Table 1. > Mechanical > Plasticity > Plastic > 
Enter Yield Stress and Plastic Strain in 
accordance with Table 1. > Click OK. 


li Edit Material x 


Name: Al-T3 


Material Behaviors 





General Mechanical Thermal Other Delete | 
Plastic 
Hardening: [Isotropic 7] nd Suboptions | 


T Use strain-rate-dependent data 


[” Use temperature-dependent data 


Number of field variables: 0 4 


Data 


Plastic 





Strain 
1 3.47E+008 it} 


2| 5.699E+008 0.1577 


Cancel 





Edit Material Property Dialog Box. 
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Section > Create > Select Shell Category > 
Homogeneous Type > Click OK. > In Edit Section 
window, enter shell thickness 0.0015875 m > 
Select Material Al1-T3 > Click OK. (Create second 
section “Section-2” with shell thickness 0.003175 
m) . 

Assign > Section > Select both stiffener joints 
in viewport > Click OK > Select “Section-2”. 
Assign “Section-1” to the rest of plate using the 
same procedure. 


5s Parts assembly is accomplished next. This is 
performed in the Assembly Module by creating an 
instance of each part in model and reorient to the 


desired position. From main menu bar: 


Instance > Create > Select Rig Ball > Click 
Apply > Select Stiffened_ Panel 








Figure 26. Ball & Stiffened Plate Impact Analysis Setup. 
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Analysis step is created next via the Step Module. 





This part allows selection of type of analysis 








needed to perform. In this case, dynamic explicit 


analysis is chosen. From the main menu bar: 


Step > Create > General procedure type > Dynamic 
Explicit > Continue > In Edit Step window, Click 
OK. 


Output can be requested in the Step Module. For 
this analysis, the default output request is used. 
To rename the both the history output and field 


output request: 


Output > Field Output Requests or History Output 
Requests Rename > Enter name or use default. 


Contact and interaction property are defined in 








the Interaction Module. From main menu bar: 


Interaction > Property > Create > Contact type 
> Continue > Mechanical Tangential Behavior > 
Frictionless > Mechanical ~> Normal Behavior > 
“Hard” Contact > Allow separation after impact > 
Click OK. 


Interaction > Create > Surface-to-surface 
contact (Explicit) > Continue > Follow direction 
at bottom of viewport > select plate as first 
surface > select rigid ball as second surface > 
Click Done. 


Boundary conditions and initial condition are 


defined in the Load Module. The plate is rigidly 





constrained on four edges. The reference point on 





the ball is also rigidly constrained except for 





the direction of motion of impact (U2). The 
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2 


ae 


Figure 27. 


10. 


initial condition or velocity of Rigid ball is 
applied at the ball’s reference point. From the 


Main Menu Bar: (Figure 27) 


BC > Name Plate Edge BC > Select Step-1 73 
Category Mechanical > Types Displacement/Rotation 
> Continue > Follow instruction at bottom of 
viewport > select all four edges > In Edit 
Boundary Condition Box, select all 6 DOFs (3 
displacements/ 3 rotations) > Click OK 


Predefined Field > Create > Initial Step 7> 
Category Mechanical > Types Velocity > Continue 
> Select reference point (RP) on Rigid Ball > 
Enter velocity for U2 (10, 30, 35, 55, 60, 65, or 
70m/s) . 








Applied Boundary Condition & Initial Velocity 








In order to mesh the fuselage panel, it must first 





be seeded. There are various ways to seed, but 


seed by “Edge Size” will be used here. Once the 





panel is seeded, the model can be selected to 
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mesh. Depending on what kind of instance was 


created in the Assembly Module (Dependent or 





Independent), the model must be meshed as such. 





If the instance was Dependent, then the model must 








be meshed via the original part because the 





instance depends on original part. If instance 





was Independent, then model must be meshed via 
assembled part because it does not depend on the 


original part. A part modeled as rigid (rigid 





ball) can not meshed in ABAQUS. In this case, we 


seed the flat panel separately from the stiffener 





support. In the main menu bar: 


Seed > Edge by size > select an edge of 
stiffener to create local seed > Click Done ~> 
Enter seed size 0.00635 m (Figure 28). 





Figure 28. 


Edge of panel seeded at 0.00635 m. 


Seed > Edge by size > Select edge of stiffener 
support to create local > Click Done > Enter 
size of seed 0.00635 m (Figure 29). 
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Figure 29. Edge of panel and stiffener support seeded at 
0.00635 m. 


Mesh > Instance > Select the whole panel 
instance to mesh > Click Done (Figure 30) 





2 Op 
: © 
oe ) 
Figure 30. Meshed stiffened panel at 0.00635 m seed. 


11. The last step is to create a job to be submitted 
for analysis using the Job Module. ABAQUS allows 
different types of job submission such as full 
analysis, restart, recover (from a crash) or 
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continue analysis. It also allows configuration 





of computer memory usage for analysis as well as 
option for multiple processors to run in parallel. 


From menu bar: 


Job > Create > Submission tab, select Full 
analysis > Memory tab, enter 512 for both > 
Precision tab > ABAQUS/Explicit precision double 
> Nodal out precision full > Click OK. (Figure 





31) 
Mm Edit Job x 
Name: StiffenedPlate_Center_Impact 
Model: Stiffened Plate_Center_Impact 
Description: [enter impact analysis of stiffened plate at 10mjs usin 
Submission | General | Memory | Parallelization | Precision | 
Job Type 
@ Full analysis ™ Recover (Explicit) 
™ Data check Restart 
™ Continue analysis 
Run Mode 
@ Background © Queue: 
Submit Time 
@ Immediately 
O wo hrs. eal min. 
{ Tip... | 
Cancel | 
Figure 31. Create a job for analysis. 


12. 


At this point an input file should’ve been 
created and can be submitted to analysis. The 
different impact analysis as listed in Table 3 can 


be performed for unstiffened panel by simply 





changing the initial velocity and resubmit the 


job. 
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APPENDIX B. STIFFENED AND UNSTIFFENED PANEL 
EIGENVALUE CALCULATION AND COMPRESSION TEST 
MODELING USING ABAQUS/CAE 


This appendix explains in detail how models’ for 


eigenvalue calculation analysis and compression test 





analysis were created using ABAQUS/CAE. The final result 

















would be an input file which is used to submit into a job 


request for analysis using ABAQUS/Standard. ABAQUS is an 





advanced finite element analysis that provides complete and 


powerful solutions for linear and nonlinear ngineering 








problems. It is a suite of finite element analysis modules. 





ABAQUS/CAE, having a modern graphical user interface (GUI) 
of menus, icons, and dialog boxes, provides the most 


complete interface with the ABAQUS solver programs available 








[Ref. 2]. In order to understand the process, the user 
should be familiar with the components of the ABAQUS/CAE and 


the appearance of the window. Figure 20 shows the components 





that appear in the main window. 


A. EIGENVALUE CALCULATION/BUCKLING PREDICTION ANALYSIS 
MODELING USING ABAQUS/CAE 


ae The part (deformed stiffened plate) involved in 
this analysis was imported from the impact 
analysis created in Appendix A. into a new model 
using the *import option. From the main menu bar 


create a new model and: 
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a 


Figure 32. 


File > Import > Part > File Filter, select 
Output Database (*.odb*) > select output file 
IMPACT StiffenedPlate C_v35.odb > Click OK > 
Select STIFFENED PANEL-1 > Rename part name 
“Stiffened_ panel” > Select “Import deformed 
configuration.” > Click OK. (Figure 32) 








Imported deformed fuselage panel at 35m/s 
impact speed. 


Once the deformed part has been imported, it has 





to be assembled in the Assembly Module. Even 
though only one part is involved and instance has 
to be created. Ensure instance name in this model 
is same as instance name in the previous model. 


From main menu bar: 


Instance > Create > in Create Instance box, name 
instance “Stiffened Panel-1 > Click OK. 
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Next, an analysis step must be created in the Step 


Module. Here, the Linear Perturbation Buckle 





procedure is used. Output of the analysis is also 











requested in this module. In the main menu bar: 


Step > Create > Select Linear Perturbation 
procedure type > Choose Buckle > Click OK > in 
Edit Step box, select Subspace Eigensolver > 
Request 5 eigenvalues > Click OK. 


Output > Field Output Request > Edit > F-Output- 
1 > Click OK to accept the default request. 


Next, the imported geometry must be defined as 








initial condition £Eor this new model. In 


addition, boundary conditions and load must also 





be defined as well. Thes ar don in the Load 





Module. From the main menu bar: 


Predefined Field > Choose Step “Initial” > 
Category “Other” > Type “Initial State” > Click 
Continue > Select the fuselage panel instance in 
the viewport to be assigned initial state > 
Bottom of viewport click Done > in Edit 
Predefined Field, enter the previous analysis job 
name “IMPACT StiffenedPlate _C_V35” > Click OK. 
(Figure 33) 
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Figure 33. 


li Edit Predefined Field x 


Name: Predefined Field-1 
Type: Initial State 










Step: Initial 
Instances: (Picked) 


Import Criteria 


Job name: [IMPACT_StiffenedPlate_C_¥35 
Step: @ Last © Specify: 
Frame: @ Last © Specify: | 


[” Update reference configuration 


OK | Cancel | 





Specify exact job name from previous impact 
analysis. 


BC > Create > Name Simple Support > Select Step- 
1 > Category “Mechanical” > Type 
“Displacement/Rotation” > Click Continue 73 
Select the far edge from screen to applied BC; may 
need to rotate and zoom in to select all nodes on 
this edge> at bottom of viewport, click Done > 
in Edit Boundary Condition box, set 
U1=U2=U3=UR2=UR3=0 (simple support) > Click OK. 
(Figure xxx) 


BC > Name “Compressed Edge” > Step-1 > Category 
“Mechanical” > Types “Displacement/Rotation” > 
Continue > Set BC for closest edge on screen 
U1=U2=UR2=UR3=0 > Click OK. (Figure 34) 
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@ Edit Boundary Condition x! 


Name: Simple Support 
Type:  Displacement/Rotation 
Step: Step-1 (Buckle) 


Region: 


Use BC for 
@ Stress perturbation and buckling mode calculation 
© Stress perturbation only 


Buckling mode calculation only 


csys: (Global) Edit... 


Distribution: 


IV ut: 
MV ue: 
MV us: 
T uri: 
i UR2: 
I¥ URS: 


Note: The 


maintained in subsequent steps. 


OK | 





Figure 34. 





@@ Edit Boundary Condition x! 


Name: Compressed Edge 




















Type: Displacement/Rotation 
Step: Step-1 (Buckle) 
Region: (Picked) 


Use BC for 
@ Stress perturbation and buckling mode calculation 
© Stress perturbation only 


© Buckling mode calculation only 


csys: (Global) Edit... 


Distribution: [Uniform = ttiti<‘isé«~S id 
I ut: jf © 
I u2: fo = =—t—~— 
Tus: ————————— 
T uri: | radians 
IV uR2: fo =83FSti<‘iCS;~;~™ radians 
I URS: fo radians 


[Uniform Ba 


| radians 
fo radians 
fo radians 


Note: The displacement value will be 
maintained in subsequent steps. 


displacement value will be 





Cancel | 


Cancel | 


Boundary condition for fixed end and loaded 
end. 


Load > Create > Step-1> Category “Mechanical” > 
Type “Shell edge load” > Continue > Select same 
edge as the “Compressed Edge” BC > Click Done at 
bottom of viewport > in Edit Load box, set 1 
Newton for “Magnitude” > Click OK. Figure (35 & 
36) 


Last step is to create a job for eigenvalue 
calculation using the Job Module as shown in 
Appendix A. For this job, use default job 


settings. 


The model can be submitted for analysis after the 


job is created. 
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Mm Edit Load x 


Name: Compressive Load 










Type: Shell edge load 
Step:  Step-1 (Buckle) 
Region: (Picked) 


Distribution: [Uniform >| Create... 
Traction: [Normal >| 
Magnitude: Fi 


Traction is defined per unit [deformed area >| 


Follow rotation 
OK | Cancel | 


Figure 35. Apply unit load for eigenvalue calculation. 








as 


Figure 36. Boundary conditions and unit load applied. 
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B. COMPRESSIVE TEST ANALYSIS MODELING USING ABAQUS/CAE 


1. Create a new model. From the main menu bar: 


Model > Create > Name 
“StiffenedPlate C_V35 PostBuckle” > Click OK. 








Bit Import the deformed geometry from impact analysis 


using the same procedure as previous model. 


3 Create a part instance in the Assembly Module 


using the same procedure as previous model. 


4. Create analysis step-l in the Step Module. Accept 


the default output request. 


Step > Create > General procedure type > Select 
Static, Riks method > Continue > Nlgeom set to 
“ON” > Click OK (Figure 37) 


Mi Create Step x 
Name: [step-1 


Insert new step after 


Procedure type: [General Ba 


Dynamic, Explicit 

Dynamic, Temp-disp, Explicit 
Geostatic 

Heat transfer 

Mass diffusion 


Soils 
Static, General 


Continue... | Cancel | 


Figure 37. Static Riks method used in analysis step. 
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Figure 38. 





In the Load Module, create boundary conditions as 
in the Eigenvalue Calculation analysis. The 


imported deformed geometry and material state is 














set aS initial similar to previous model. Initial 





load used during Riks method is 50000 Newton. In 


the main menu bar: 


Load > Name “Compressive Load” > Step-1 73 
Category “Mechanical” > Type “Shell edge load” > 
Continue > Select closest edge in viewport > 
Click Done > in Edit Load box, enter 50000 Newton 
for “Magnigtude” > Click OK. (Figure 38) 


lm Edit Load x 


Name: Compressive Load 
Type: Shell edge load 
Step:  Step-1 (Static, Riks) 
Region: (Picked) 


Distribution: [Uniform 7] Create... | 
Traction: [Normal 7] 
Magnitude: [s0000 


Traction is defined per unit [deformed area 7] 


Follow rotation 





_ «| 


Apply load for Riks buckling method. 


Create job for analysis similar to previous model. 





Use default setting for this job. 
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Once the job has been created, It can be submitted 





for a full analysis. All compressive failure 
models can be created in a similar fashion by 


following the above listed steps. 
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